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The Effects of Galaxy Shape and Rotation on the X-ray Haloes of 
Early-Type Galaxies 
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ABSTRACT 

We present a detailed diagnostic study of the observed temperatures of the hot X-ray coro- 
nae of early-type galaxies. Extending the investigation carried out in Pellegrinil (1201 ll) with 
spherical models, we focus on the dependence of the energy budget and temperature of the 
hot gas on the galaxy structure and internal stellar kinematics. By solving the Jeans equations 
we construct realistic axisymmetric three-component galaxy models (stars, dark matter halo, 
central black hole) with different degrees of flattening and rotational support. The kinematical 
fields are projected along different lines of sight, and the aperture velocity dispersion is com- 
puted within a fraction of the circularized effective radius. The model parameters are chosen 
so that the models resemble real ETGs and lie on the Faber-Jackson and Size-Luminosity re- 
lations. For these models we compute T„ (the stellar heating contribution to the gas injection 
temperature) and T„ (the temperature equivalent of the energy required for the gas escape). 
In particular, different degrees of thermalisation of the ordered rotational field of the galaxy 
are considered. We find that T„ and r„ can vary only mildly due to a pure change of shape. 
Galaxy rotation instead, when not thermalised, can lead to a large decrease of T*; this effect 
can be larger in flatter galaxies that can be more rotationally supported. Recent temperature 
measurements Tx, obtained with Chandra, are larger than, but close to, the r« values of the 
models, and show a possible trend for a lower Tx in flatter and more rotationally supported 
galaxies; this trend can be explained by the lack of thermalisation of the whole stellar kinetic 
energy. Flat and rotating galaxies also show lower Lx values, and then a lower gas content, 
but this is unlikely to be due to the small variation of T„ found here for them. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: fundamental parameters - galax- 
ies: ISM - galaxies: kinematics and dynamics - X-rays: galaxies - X-rays: ISM 
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1 INTRODUCTION 

High-quality X-ray obsei^ations of early-type galaxies (ETGs) per- 
formed with the Chandra X-ray Observatory have produced a large 
body of data for the study of the hot gas haloes in these galaxies 
with unprecendented detail. In particular, the nuclear and the stellar 
(resolved and unresolved) contributions to the total X-ray emission 
could be subtracted, obtaining more accurate properties of the hot 
interstellar medium (ISM) than ever before. From a homogeneous 
and thorough X-ray analysis for the pure gaseous component, sam- 
ples of ETGs have been built with an improved measurement of 
the X-ray average temperature Tx and luminosity Lx for the gas 
only (e.g.. ISoroson et al.11201 ll hereafter BKF). This X-ray infor- 
mation can now be compared with that provided by other funda- 
mental properties of ETGs, such as their total optical luminosity 
(Lb or Lk in the B or K band), central stellar velocity dispersion 
{(T^), galaxy rotation and shape, to explore and revisit the possible 
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relations betweeen the main galactic properties and those of the X- 
ray gas. For example, in the Lx - Lk correlation, the previously 
known large variation of up to two orders of magnitude in Lx at the 
same Lk has been even extended, due to the inclusion of hot-gas 
poor ETGs in a larger fraction than previously possible (with Lx 
extending down to ~ 10^* erg s"'; BKF). This variation is related 
to the ISM evolution over cosmological time-scales, during which 
stellar mass losses and explosions of Type la supemovae (SNIa) 
provide gas and gas heating, respectively, to the hot haloes (possi- 
bly in conjunction with feedback from accretion on to the central 
supermassive black hole, SMBH). Modulo environmental effects 
like galaxy interactions or tidal stripping, the hot gas content and 
temperature fundamentally depend on the energy budget of the hot 
ISM, that in turn depends on the particular hos t galaxy structure 
and internal kin ematics (e.g.. Iciotti et al.lll99lL hereafter CDPR; 
|Pellegrinil201ll . hereafter PI 1). For example, the luminous and dark 
matter content and distribution determine the potential well shape 
and depth, and so the binding energ y of the gas and its dynamical 
state (e.g.. lCiotti & Pellegrini 1 19961 hereafter CP96). Indeed, one 
of the discoveries that followed the analysis of first X-ray data of 
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ETGs was the sensitivity of tiie hot gas content to major galaxy 
properties as the shape of the mass distributio n, and the mean rota- 
tion velocity of the stars (see |Pellegrinill20 1 21 for a review). The in- 
vestigation of the origin of this sensitivity is the goal of the present 
paper. 

A relation between the hot gas retention capability and the 
intrinsic galactic shape became apparent already in the X-ray sam- 
ple of ETGs built from Einstein observations: on average, at any 
fixed Lb, rounder systems had larger total Lx and Lx/L^, a mea- 
sure of the galactic hot gas content, than flatter ETGs and SO 
galaxies JEskridge et al.lll995h . Moreover, galaxies with axial ratio 
close to unity spanned the full range of Lx, while flat systems had 
Lx S 10*' erg s"'. This result was not produced by flat galaxies hav- 
ing a lower Lb, with respect to round ETG s, since it held ev en in the 
range of Lb where the two shapes coexist jPellegrinil 19991) . The re- 
lationship between Lx and shape was reconsidered, confirming the 
above trends, for the ROSAT PSPC sample (Pellegrini 2012), and 
for the Chandra sample JLi et al.l2011ah . Therefore, there seems to 
be an empirical dependence of the hot gas content on the galactic 
shape, and it was suggested that a flatter shape by itself may be 
linked to a less negative binding energy for the gas (CP96). How- 
ever, since flatter systems also possess a hi gher rotational support 
on average (e.g., Isirmev & Tremainelll987h . also the influence of 
galactic rotation on the hot gas was called into question. For ex- 
ample, in rotationally supported ETGs the gas may be less bound, 
compared to the ISM in non-rotating ETGs, leading rotating ETGs 
to be more prone to host outflowing regions. For these reasons, the 
efl^ects on Lx of both galactic shape and rotation were studied for 
a sample of 52 ETGs with known maximum rotational velocity of 
the stars Vmax, and so wit h a measure of Knax/o "c, an indicator of 
the importance of rotation (Pellegr ini et al.ll997h . It was found that 
Lx/Lb can be high only for Vmax/o"c < 0.4, and is limited to low 
values for Vmax/Cc > 0.4. This trend was not produced by being 
the ETGs with high Vniax/cTc confined to low Ln. lSarzi et alj J2010h 
investigated again the relationship between X-ray emission (from 
Einstein and ROSAT data) and rotational properties for the ETGs 
of the SAURON sample, confirming that slowly rotating galaxies 
can exhibit much larger luminosities than fast-rotating ones. 

Renewed interest in the subject has come recently after the 
higher quality Chandra measurements of Lx and Tx have become 
available. In an invest igation using Cha ndra and ROSAT data for 
the ATLAS^° sample. ISarzi et alj J2013L hereafter S13) found that 
slow rotators generally have the largest Lx and Lx/Lk values, and 
Tx values consistent just with the thermalisation of the stellar ki- 
netic energy, estimated from cTe (the stellar velocity dispersion av- 
eraged within the optical eff^ective radius R^]. Fast rotators, instead, 
have generally lower Lx and Lx/Lk values, and the more so the 
larger their degree of rotational support; the Tx values of fast rota- 
tors keep below 0.4 keV and do not scale with cr^ (see also BKF). 
Considering that fast rotators are likely to be intrinsically flatter 
than slow rotators, and that the few slow rotators with low Lx are 
also relatively flat, S13 supported the hypothesis whereby flatter 
galaxies have a harder time in retaining their hot gas (CP96). To 
explain why fast rotators seem confined to lower Tx than slow rota- 
tors, they suggest that the kinetic energy associated with the stellar 
ordered motions may be thermalised less efficiently. 

In order to help clarify what is the expected variation of hot 
gas content and temperature, originating in a variation of shape and 
internal kinematics (and in its degree of thermalisation) of the host 
galaxy, we embarked on an investigation based on the numerical 
building of state-of-the-art galaxy models, and of the associated 
temperature and energy budget for the gas. In Sect.|2]we define a set 



of mean temperatures for the models, some of which already intro- 
duced in PlI. In Sect.[3]we describe the different profiles adopted 
for the mass components of the galaxy models, the scaling laws 
considered to constrain the models to resemble real galaxies, the 
observable properties of the models in the optical band, and the 
procedure to obtain flat models. Our main results are presented in 
Sect.|4] together with a comparison with the observed X-ray prop- 
erties of ETGs. Finally Sect. |5]presents our main conclusions. Ap- 
pendix|A]recalls the fluid equations in the presence of source terms, 
and AppendixiBlsummarizes the main procedural steps of the code 
built on purpose for the construction of the models. 



2 THE TEMPERATURES 

Here we introduce a set of gas mass-weighted temperatures, equiv- 
alent to the injection and binding energies of the hot gas in ETGs. 



2.1 The injection temperature 

In the typically evolved stellar population of ETGs, the main pro- 
cesses responsible for the injection of gas mass, momentum and en- 
ergy in the ISM are stellar winds from red/asymptotic giant branch 
stars, and Type la supernova explosio ns (SNIa), the only one s ob- 
served in an old stellar population (e.g.'Cappel laro et al.ll999h . The 
wind material outflowing from stars leaves the stellar surface with 
low temperatures and low average velocities (~ few 10 kms"'; 
IParriott & BregmaiJuOOSh , so that all its energy essentially comes 
from the stellar motion inside the galaxy. SNIa's explosions, in- 
stead, provide mass to the ISM through their very high velocity 
ejecta (~ few IC* km s"'). Thus, this material provides mass and 
heat to the hot haloes, via thermalisation of its energy through 
shocks with the ambient medium or with other ejecta, heating up 
to X-ray emitting temperatures. 

The injection energy per unit mass due to both heating pro- 
cesses is ginj = 3^B7'mj/(2yUOTp), whcre ^b is the Boltzmann constant, 
fi = 0.62 is the mean molecular weight for solar abundance, mp is 
the proton mass, and Linj is defined as 



M,T, + MsnLsn 
M 



(1) 



Here T, and Lsn are the injecta temperatures resulting from the 
thermalisation of their interactions with the ISM through stellar 
winds and SNIa's respectively (see below). M is the total mass- 
loss rate for the entire galaxy, given by the sum of the stellar 
mass-loss rate M, and of the rate of mass loss via SNIa events 
MsN (^ = M, + Msn). The time evolution of the stellar mass- 
loss rate M, can be calculated using single-burst stellar popula- 
tion synthesis models for different init ial mass functions (IMFs) 
and metallicities (e.g.. iMarastonI uOOSh . For example, at an age 
of 12 Gyr, M.iM^ yr ~') ~ 2x I0~" Lb(Lbo) for the Salpeter 
or Kroupa IMF (e.g.. |Pellegrinill20I2l) . Msn is instead given by 
MsN = A^SN^SNi where Msn = 1-4 Mq is the mass ejected by 
one SNIa event and i?sN is the explosion rate. For local ETGs it is 
RsN = 0.16(//o/70)- X 10~'^ LB(LB 0)yr~', where Hpis the Hubble 
constant in units of km s^'Mpc"' JCappellaro et al.lll999 V More 
recent measure ments of the ob served rates of supemovae in the 
local universe ( ILi et al. 201 Ibl) give a SNIa rate in ETGs consis- 
tent with that of lCappellaro et al.1 ( Il999l) . For this rate, and Hq =70 
kms^'Mpc"', one obtains Msn = 2.2 x 10"'^Lb(Lbo) Mgyr"', 
which is ~ 80 times smaller than the M, above for an age of 12 
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Gyr. Thus the main source of mass is provided by M,, and approx- 
imating M ^ M,, we have Tinj ^T, + {M^^IM,)T^^. 

Neglecting the internal energy and the stellar wind velocity 
relative to the star, T, is the sum of two contributions, deriving 
from the random and the ordered stellar motions. In axisymmetric 
model galaxies as built here, the stellar component of the galaxy 
is allowed to have a rotational support, and the latter can be con- 
verted into heating of the injected gas in a variable amount. The 
extent of the contribution of rotational motions is not known a pri- 
ori, since it depends on both the importance of the stellar ordered 
motions and the dynamical statu s of the surroun ding gas al ready 
in situ (see Appendix |AJ see also D'Ercole et al.l l2000: Negri et alj 
l2013ah . Given the complexity of the problem, hydrodynamical sim- 
ulations are needed to properly calculate this heating term, but we 
can still obtain a simple estimate of it by making reasonable as- 
sumptions. We define the equivalent temperature of stellar motions 
T, as 



To- + TthT'rot 



where 



fimp 



- Jp,Tr(a' 



SknM. 
is the contribution of stellar random motions 



)dV 



T,. 



3kjiM. 



/p.v/ 



dV 



(2) 



(3) 



(4) 



is the one due to the stellar streaming motions, and yth is a parame- 
ter that regulates the degree of thermalisation of the ordered stellar 
motions. M, is the stellar mass of the galaxy, cr^ is the velocity 
dispersion tensor, and v^ is the azimuthal component and the only 
non-zero component of the streaming velocity v = v^e^ (see Ap- 
pendix|Bj. In Eqs. JSJ and lO, as in the remainder of the paper, we 
assume that the gas is shed by stars with a spatial dependence that 
follows that of the stellar distribution p., so that the density pro- 
file of the gas injected per unit time is proportional to p., i.e. it is 
^ = MpJM, in Eqs. l lA2t - l |A3t of Appendix lAl 
The parameter 7,1, is defined as 



7ih 



pmp 



r,.o,3A:BM, 



/' 



ii«-virdK 



(5) 



where u is the velocity of the pre-existing gas (see Appendix|Aj. A 
simple estimate for 7,1, is obtained when the gas velocity is propor- 
tional to V, i.e., u = av, where a is some constant. In this special 
case, from Eqs. (|4) and (|5} it follows that 7th = (or - 1)^. When 
a = I, gas and stars rotate with the same velocity and no ordered 
stellar kinetic energy is thermalised, whereas for a = the gas is 
at rest and all the kinetic energy of the stars, including the whole of 
the rotational motions, is thermalised. Clearly both cases are quite 
extreme and unlikely, and plausibly the pre-existing gas will have 
a rotational velocity ranging from zero to the streaming velocity of 
stars, i.e., < a < 1 and then 1 > yth > (the case of a constant 
a > 1, where the pre-existing gas is everywhere rotating faster than 
the newly injected gas is not considered). Note that, contrary to 
TthT'rot, r^ is in principle exact, and can be computed a priori. 

The internal plus kinetic energy of the ejecta, released during 
a SNIa event, is of the order of iisN - 10" erg. Depending on the 
conditions of the environment in which the explosion occurs, the 
radiative losses from the expanding supernova remnant may be im- 
portant, and so the amount of energy transferred to the ISM through 
shock heating is a fraction rj of iis^. Realistic v alues of rj for the hot 
and diluted ISM of ETGs are around 0.85 (e.g.. lTang & Wan J2005l 



lThomtonet"ai]|l998h . thus 



TsN — 



2yump TjEs 
3^B Msi 



(6) 



and, substituting the above expressions for Msn and M,, we obtain 
the average injection temperature 

1 



T, + 1.7—^ X 10' K. 
0.85 



(7) 



A possible additional source of heating for the gas could be 
provided by a central SMBH. Through its gravitational influence, it 
is responsible for the increase of the stellar motions within its radius 
of infl uence (of the order of a few tens of parsecs; e.g., Pelle grinil 
12013)). We consider this effect here, while we neglect possible 
effects as radiative or mechanical feedback. 



2.2 The temperatures related to the potential well 

The gas ejected by stars can be also heated 'gravitationally' by 
falling into the galactic potential well to the detriment of its po- 
tential energy, and by the associated adiabatic compression. When 
stellar mass losses accumulate, the gas density can reach high val- 
ues and the cooling time can become smaller than the galactic age; 
if the radiative losses increase considerably, the gravitational force 
overwhelms the pressure gradient and eventually the gas starts in- 
flowing toward the centre of the galaxy. Thus we can define a tem- 
perature 

2fimpE+ 2pmp 



t: 



■I 



p. (1'-*o)dV 



(8) 



3ytB 3ksM, 

where E'^ is the average change in gravitational energy per unit 
mass of the gas flowing in through the galactic potential 0(x) down 
to the galactic centre, and Oq = 'l'(O). Note that E* > 0, having 
assumed as usual that <I>(x) < 0. However, as discussed in Pll, 
most of £^ may be radiated away, and there are conditions under 
which Eq. l[8j does not apply. Therefore, given these uncertainties, 
we consider T* just as a reference value, and keep in mind that 
the temperature achievable from infall can be much lower than that 
given by Eq. itSj. 

By analogy with r+, we can define a temperature 



2fimpE„ 

3kn 



Ifinif 



fp. <I>dV, 



(9) 



3A:bM, 

where E^ is the average energy necessary to extract a unit of gas 
mass from the galaxy, with the assumption that 0(00) = 0. If the 
gas rotates, Eq. (|9) must be modified since, thanks to the centrifugal 
support, the gas is less bound. Assuming again that u = av, then 



E,(a) : 



-wj' 



^ « — - 
O + y V/ 



dV, 



(10) 



so that T„(a) = T^ - a^Ty^f Note that, for a given r^o,, the smallest 
is 7,h (the largest is a), the smallest is T, (the gas is less heated), 
but also the lower is T^ (the gas is less bound). 

When the galaxy mass distribution has a potential that di- 
verges at small and/or large radii, as for the singular isothermal 
sphere, we assume the gas has been extracted from the galaxy when 
it has reached a distance of 15 i?e from the galactic centre, so that 
E'^ and E'^ do not correspond exactly to Eqs. dlO) and ([8}. Finally, if 
energy losses due to cooling are present, the gas would need more 
than E„ to escape, but these losses are negligible for outflows that 
typically have a low density. 

In case of gas escape, we can introduce another mass- 
weighted temperature by considering the enthalpy per unit mass 
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of a perfect gas h = yk^T/ljimpiy - 1)] = cl/(y - 1), where y 
is the ratio of the specific heats and c, is the sound speed. Build- 
ing on the Bernoulli theorem, we can derive a fiducial upper limit 
to the temperature of outflowing gas. For a fixed galactic potential, 
the energy of the escaping gas can be divided between kinetic and 
thermal energy with different combinations. In the extreme case in 
which the gas reaches infinity with a null velocity and enthalpy, and 
it is injected with a (subsonic) velocity u = av, from the Bernoulli 
equation h(x) + v^(x)/2 + 0(x) = 0, we derive a characteristic gas- 
mass averaged escape temperature 



2yump 
'SkeM, 



|p.(o+^T 



dV = -T^(al 



(11) 



for a monoatomic gas (see Pll for more details). In the opposite 
case of an important kinetic energy of the flow, the gas temperature 
will be lower than T^"l^. 

In summary, T^ is a temperature equivalent to the energy re- 
quired to extract the gas, while T™^!" is close to the temperature we 
expect to observe for outflowing gas. For inflowing gas, we expect 
to observe a temperature much lower than T^, since more than 
~ 0.5E^ is radiated away or goes into kin etic energy of the gas , 
or because of condensations in the gas (e.g.. lSarazin & Ashelll989f) . 
For reference, for realistic spherical models, E'!; ~ 2E^ , thus the 
temperature of the inflowing gas should be lower than T„ (Pll). 

Finally, we mention about the relation between observed tem- 
peratures Tx and the average mass-weighted temperatures of this 
Section. The latter are derived under the assumption that the gas 
density p„.,^ follows that of the stars, which is appropriate for the 
continuously injected gas (e.g., for T, and Tinj), while the bulk of 
the hot ISM may have a different distribution; thus, mass-weighted 
T„ and T^^* referring to the whole hot gas content of an ETG 
may be different from those given by Eqs. (|9) an d (lilt. In general, 
the Pgas profile i s shallower than that of p, (e.g.. lSarazin & White! 
1 19881 : iFabbiandl 1989.) . and then the gas mass-weighted T^ would 
be larger than derived with Eq. (|8), and the mass-weighted T~ or 
r™^!' would be lower than derived using Eqs. (|9) and (lilt . For steady 
winds, instead, when the gas is continuously injected by stars and 
expelled from the galaxy, the Pg^s oc p, assumption is a good ap- 
proximation. 

Another point is that the Tx values are emission-weighted 
averages, and will coincide with mass-weight ed averages only if 
the entire ISM ha s one temperature value (e.g.. lCiotti & Pellegrinil 
120081 : lKimll20I2h . A single Tx value measured from the spectrum 
of the integrated emission will tend to be closer to the tempera- 
ture of the densest region, in general the central one, thus it will 
be closer to the central temperature than the mass-weighted one. 
The temperature profiles observed with Chandra tend to be quite 
flat, except for cases where they increase outside of ~ 0.5i?e (gen- 
erally in ETGs with the largest Tx), and for cases of negative tem- 
perature gradients (in ETGs with the lowest Tx', Diehl & Statleg 
l2008l : iNagino & Matsushitall2009l) . Therefore, the largest Tx may 
be lower than mass-weighted averages, and the lowest Tx may be 
larger than them. In conclusion, the comparison of Tx and the gas 
content with the gas temperature and binding energy introduced in 
this Section (as T, and E„) represents the easiest approach for a gen- 
eral, systematic investigation involving a wide set of galaxy models, 
but the warnings above should be kept in mind. Note, however, that 
the conclusions below remain valid when taking into account the 
above considerations. 



3 THE MODELS 

The galaxy models used for the energetic estimates include three 
components: a stellar distribution, a dark matter (DM) halo, and 
a central SMBH. The stellar component is axisymmetric and can 
have different degrees of flattening, while for simplicity the DM 
halo is spherical. The SMBH is a ce ntral mass concentration with 
mass Mbh = lO^^'M., following the iMagorrian et alj ( Il998h rela- 
tion. Its effects are minor, but it is considered for completeness. 
For these models, the Jeans equations are solved under the stan- 
dard assumption of a two-integral phase space distribution function 
(see Appendix IB]). so that, besides random motions, stars can have 
ordered motions only in the azimuthal direction. The decomposi- 
tion of the azimuthal motions in velocity dispersion and stream- 
ing velocity is performed via the A^-decomposition introduced by 
ISato h ( 1980): thus, the amount of rotational support is varied sim- 
ply through the parameter k. With the adoption of the mass profiles 
detailed below for the stars and the DM, we built galaxy models that 
reproduce with a good level of accuracy the typical properties of 
the majority of ETGs. The models are then projected along two ex- 
treme lines of sight (corresponding to the face and edge-on views) 
and forced to resemble real galaxies as described in Section [331 



3.1 Stellar distribution 

The stellar distribution is described by the de Vaucouleurs] ( 119481) 
law, by using the deprojection of iMellier & Mathez ( 1987i) general- 
ized for ellipsoidal axisymmetric distributions 



MR,z)=por°'''^^P(-e'\ 



with 



M,b' 



Pa ■ 



l6nqRl^r(8.5S)' 



IM- 



(12) 



(13) 



where (R, if, z) are the cylindrical coordinates and b ^ 7.66925. The 
flattening is controlled by the parameter g < 1, so that the minor 
axis is aligned with the z axis. R^o is the projected half mass radius 
(effective radius) when the galaxy is seen face-on: for an edge-on 
view, the circularized effective radius is R^ = R^-g yfq (Sect. 13.41 
and Appendix iBt. We assume a constant stellar mass-to-light ratio 
T, all over the galaxy, so that M, is directly proportional to the 
luminosity L. Note that Eq. ( I13t guarantees that the total stellar 
mass (luminosity) of the model is independent of the choice of q 
and/feo- 



3.2 Dark matter halo 

Given the uncertainties affecting our knowledge of the density pro- 
file of DM haloes, we explored four families of DM profiles. The 
first one is the scale-free singular isothermal sphere (SIS) 



Ph(r) 



4nGr'- ' 



Oh(r) = v-lnr. 



(14) 



where v^ is the halo circular velocity. The gravitational potential of 
this profile diverges at small and large radii, thus the potential is 
truncated at a distance of 15 ^e to obtain a finite T^ 



A number of recent works are reconsidering the lEinastol 

1 1965h profil e as ap propriate to model DM haloe s (e.s 



Navarro et al. 



as ap propriate to modei DM hatoe s (e.g . 
12004 iMerritt et al.l l2006t iGao et al.l l2008l 
I2OI0I) . The density distribution of this profile 
is the three-dimensional analogue of the Sersic law, widely used to 



Navarro et al. 
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fit the surface brightness profiles of ETGs. The density is described 
by 



Ph(r) = Pc exp(rf„ - x), 



(15) 



where p^ is the density at the volume half- mass radius ?■[,, x = 
d„(r/ri,y'", n is a free parameter, and d„ is well approximated by 
the relation 



3n + , 

3 1215 n 



(16) 



jRetana-Montenegro et al.l2012n . Finally the gravitational potential 
is 

GMh 



Oh(-<:) 



r(3n,x) x'T(2n,x) 



r(3n) r(3n) 

The third family is based on the lHemquistI ( 119901) profile 



P\ir) 



2nr{r + rn)^ ' 



l'h(r) : 



r + rh 



(17) 



(18) 



where Mh and rn are the halo total mass and scale radius, respec- 
tively. 

Lastly, we used also the NFW profile JNavarro et al.lll997h 



Phir) 



Pent S^r^ 



(19) 



r(l +r/rh)- 

wherepcrit = 3H-/SnG is the critical density for closure. The total 
mass diverges, so it is common use to identify the characteristic 
mass of the model Mh with the mass enclosed within r2Qo, defined 
as the radius of a sphere of mean interior density 200 pcrit- Then, 
from the definition of ^200, the concentration c = rimh'h and the 
coefficient 6r are linked as 



6c 



200 



3 ln(l+c)-c/(l+c) 
The gravitational potential of the NFW profile is 



^h{r) = -4nGpcnt S^i^ 



Ml + r/r^,) 



(20) 



(21) 



3.3 Linking tlie models to real ETGs 

One of the most delicate steps of the present study is to have a sam- 
ple of galaxy models, characterised by various degrees of flattening 
and rotational support, that closely resemble real ETGs, at least in a 
statistical sense. This is accomplished by flattening spherical mod- 
els, that we call 'progenitors'. 

In fact, the process of flattening a galaxy model is not triv- 
ial, and it is highly degenerate, as illustrated by the exploratory 
work of CP96, where the full parameter space of two-component 
Miyamoto-Nagai models was explored. Here, we begin with a 
generic spherical galaxy model, and we impose that its effec- 
tive radius R^ and aperture luminosity-weighted velocity disper- 
sion within i?e/8, cr^g, satisfy the most important observed scaling 
laws (SLs) of ETGs, the Faber-Jackson and the Size-Luminosity 
relations. In particular, we use the Faber-Jackson and the Size- 
Luminosity relations derived in the r band for a 80000 ETGs 
drawn from Data Release 4 (DR4 ) of the Sloan Digital Sky Sur- 
vey (SDSS: Desroches et alJuOOTT) . These relations are quadratic 
best-fitting curves, with a slope varying with luminosity L,.: 



logo-ei 



-1.79-1- 0.674 log L, - 0.0234 log" L„ 



(22) 



logi^e = 1.50 -0.802 log L, -1-0.0805 log^L,, (23) 
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where cr^g and R^ are in units of km s~' and kpc respec tively, and 
L, is calibrated to the AB system JDesroches et al.ll2007n . 

In practice, we fix a value for cr^g in the range 150 km s"' < 
(T^g < 300 km s"', and then we derive L, and R^ from Eqs. ([22) 
and ( 123b . After conversion of L^ to the V-banqJ (Ly), we derive 
M, adopting a (luminosity dependent) V-band mass-to-light ratio 
Y, a ppropriate for a 12 Gyr old stellar population w ith a Kroupa 
IMF ( Maraston 2005). Following empirical evidences dBender et al.l 
1992 : Cappellari et al. 2006) . we assume that T, oc L^^*, obtaining 
3.3 < T. < 4.7. With this choice, the models need a DM halo 
to reproduce the assigned iTcg. We consider the four different fam- 
ilies of (spherical) DM haloes in Sect. 13.21 whose parameters are 
fixed to reproduce the assigned iXeg. The simplest family is that 
with the SIS halo in Eq. ( I14t . where we fix v^ so that the progenitor 
has the given cr^g. For the Einasto DM haloes, we fix n = 6 and 
Th - TR^ in Eq. ( |15t , in order to obtain r^ values in the accepted 
range for ETGs (see, e.g., Merritt et al. 2006; Navarro et al. 20100, 
and to keep low the DM fraction /dm in the central regions of the 
model (see below). Mh, the only remaining free parameter, is then 
determined by the chosen cr^g . This procedure gives Mh values that 
are ^ 20 percent larger than in the SIS case, due to the shallower 
density slope of the Einasto DM halo at small radii, which trans- 
lates into a weaker effect on the stellar random motions, and then 
into a larger DM amount required to raise the central stellar velocity 
dispersion profile up to the chosen cr^g. Also for the Hemquist and 
NFW families we choose r\, ^ TR^, and Mh is fixed to reproduce 
the assigned cr^g. For the Hemquist family, this request results in 
l.SxlO'^Mo < Mh < 4x lO'^Ma, while fo r the NFW haloes we find 
12 < c < 25 dBinnev & Tremaind [19871 : iNapolitano et alj|2009h . 
corresponding to IO'^Mq > Mh > 7.2 x IO'^Mq. For all models, 
the resulting Mh/M. ratios agree with those given by cosmological 
simulations and galaxy mass functions (Naravanan & Dave 20121) . 
A summary of the properties of some spherical progenitors, for SIS 
and Einasto DM haloes, is given in Table[T] 

An important quantity characterizing the models is the ef- 
fective DM fraction, defined as the ratio of the DM mass to 
the total mass contained within a sphere of radius R^, /dm = 
Mh(,Rc)/Mtot(Rc)- We compute /dm a posteriori, to check that it 
agrees with the values found for well studied ETGs from stel- 
lar dynamics and gravitation al lensing studies (that is, /dm ~ 
0.3; Cappellari et al. 2006i. lGerhard et al.l200lLlThomas et al.l2005l , 
iTreu & Koopmans.,2004.) . In particular, for the NFW families we 
found quite high /dm values (of the order of ~ 0.66 for the progen- 
itor) due to the larger Mh values. 



3.4 From spherical to flat ETGs 

In principle, realistic flat and rotating galaxy models could be con- 
structed with a Monte-Carlo approach, where all the model param- 
eters are randomly extracted from large ranges, the resulting mod- 
els are projected and observed at random orientations, and then 
checked against the obs erved SLs, retaining o nly those in accor- 
dance with observations ( ILanzoni & Ciottil2003i) . This approach is 
unfeasible here, becaus e the model construction is based on numer- 
ical integration (while i Lanzoni & Ciottill2003l used the fully ana- 
lytical but quite unrealistic Ferrers models), and the computational 



' The V-band luminosity Ly is computed using the standard trans- 
formation equations between SDSS magnitudes and other systems 
(http://www.sdss3.org/dr9/algorithms/sdssUBVRITransform.php 1, also as- 
suming B ~ V = 0.9 as appropriate for ETGs dPonas et al.l2007l) . 
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Table 1. Fundamental galaxy parameters for the progenitors. 
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Notes: (1) Stellar velocity dispersion, as the luminosity-weighted average within an aperture of radius Rc/8. (2) and (3): luminosities in the V band (derived 
as de scribed in Sect. 13.31 and K band, from Lk = 4Lv as appropriate for a 12 Gyr old stellar population with a Kroupa IMF and solar metallicity (Marastgij 
I2OO3). (4) Effective radius. (5) Stellar mass. (6) V band stellar mass-to-hght ratio. (7) - (9) Circular velocity, dark matter mass within a sphere of radius 15 R^ 
and dark matter fraction within R^, for the SIS halo. (10) - (1 1) Dark matter mass and dark matter fraction within Re for the Einasto halo. 



time of a Monte Carlo exploration of the parameter space would be 
prohibitively large. So we solved the problem as follows. 

We flatten each spherical progenitor, acting on the axial ratio q 
and on the scale-lenght R^o of the stellar density in Eq. J13t . while 
keeping Li, T, (and then M, and Mbh), and the DM halo the same. 
For given q and R^o, the circularized effective radius R^ depends 
on the line-of-sight (l.o.s.) direction, ranging from R^o (when the 
model is observed face-on, hereafter FO) to -s/qR^o (in the edge- 
on case, EO). Thus, a request for a realistic model is that L, and 
Re remain consistent with the observed Size-Luminosity relation, 
independently of the l.o.s. direction. In turn, also cr^s will change 
due to the flattening, both as a consequence of the choice of q and 
Rco, and of the l.o.s. inclination. 

To include all possible inclination effects, from each spherical 
progenitor, we build two sub-families of flat descendants, the FO- 
built ones and the EO-built ones. In the first sub-family, R^ is the 
same of the spherical progenitor when the flat model is seen FO; 
in the other, R^ is the same of the spherical progenitor when the 
flat model is seen EO. This implies that R^o may vary: with the 
decrease of q, ifeo remains equal to R^ of the spherical progenitor 
in the FO-built case, while R^o increases as R^/ yfq in the EO-built 
case. Therefore, in this latter case, there is a consequent expansion 
(or size increase) of the galaxy, and a decrease of the galaxy scale 
density po oc yjq in Eq. | |13I ). On the contrary, in the FO-built sub- 
family, there is a density increase as po °^ 1 \ as the galaxy is 
compressed along the z-axis. 

Then, we compute o-^g according to the procedure described 
in Appendix |B] for the range spanned by the Satoh parameter 
< fe < 1. Since the FO and EO-built models are characterised by 
different structural and dynamical properties, we must check that, 
once observed along arbitrary inclinations, the models are still con- 
sistent with the observed SLs. For example, a FO-built model, when 
observed EO, will have an R^ smaller than the progenitor, while an 
EO-built model will have a larger R^ when observed FO.Therefore, 
only galaxy models that, observed along the two extreme l.o.s. di- 
rections (FO and EO), lie within the observed scatter of R^ and (Ttg 
at fixed Lf should be retained in our study. Remarkably, all the mod- 
els constructed with our procedure have been found acceptable. 

The effects of flattening on (Teg deserve some comments. For 
the EO view, cr^g of the EO-built models decreases for increas- 
ing flattening, due to the associated model expansion; cr^g further 
decreases at increasing k, as more galaxy flattening is supported 
by ordered rotation. Also for the FO view, cr^g of EO-built sys- 
tems decreases, but independently of k (affecting only cr^, while 
(Tr = (T:^). In the FO-built models, one would naively expect an 
increase of cr^g due to the density increase (and so to the gravita- 
tional potential deepening), but this is not the case: even though 



less severely than for EO-built models, cr^g still decreases (both for 
the FO and EO views). The simplest way to explain this behaviour 
is to c onsider the FO flattening o f the fully analytical Ferrers ellip- 
soids JBinnev & Tremainelll987r) . As the flattening increases, the 
density raises, the gravitational potential well deepens, and the ver- 
tical force increases, but again the velocity dispersion drops. The 
p hysical reason, behind th e mathematics [e.g., see eqs. (C4-C11) 
in lLanzoni & Ciottil ( l2003h l. is that stars need less vertical velocity 
dispersion in order to support the decreased z-axis scale-lenght, and 
so the FO view cr^g decreases. iXeg decreases less when observed 
EO, because of the decrease of R^, that causes cr^g to be computed 
within a smaller area around the galactic centre. Of course, if the 
galaxy is not fully velocity dispersion supported, the decrease of 
(Te8 for an EO view can be even larger than for the FO one, since 
part of the stellar kinetic energy is stored in ordered motions that 
do not contribute to (Ti-g. 

These general results, obtained for realistic models, about the 
variation of cr^g in flat and rotating galaxies of fixed stellar mass, 
show that some caution should be exercised when using simple dy- 
namical mass estimators based on the velocity dispersion measured 
in the central regions of galaxies. This point is particularly relevant 
for studies of the hot haloes properties, that notoriously mainly de- 
pend on the galaxy mass (CDPR, S13). 

As a further test o f the models, we also c alculated the param- 
eter Afi, introduced by lEmsellem et al.l ( 120071) and related to the 
mean amount of stellar rotational support. Our Ajt radial profiles, 
even for the k = 1 case, are in good agre ement with the profile s 
in the ATLAS^° sample of ETGs (fig. 5 in'E msellemet alj|201lh . 
for each galaxy ellipticity. Finally, our method of definition of the 
DM halo implies a constant Mi,/M, ratio within each family, but 
not a constant /dm, that depends on q (through the variation that q 
imposes to p.), as one can see in Fig.[T]for the SIS and the Einasto 
families. Note that /dm can decrease or increase with q depending 
on the construction mode: the increase of the stellar density in the 
FO-built models results in lower /dm, since the DM halo is fixed; 
the reverse is true for the EO-built models. Moreover, the Einasto 
models have always higher DM fractions than the corresponding 
SIS ones, due to the steepness of the SIS profile at small radii, re- 
quiring less DM to reproduce the chosen value of (Teg. 



4 RESULTS 

Having built a large set of realistic galaxy models, consistent with 
the observed SLs and with DM haloes in agreement with current ex- 
pectations, we can now study the effects of flattening and rotational 
support on the temperatures of the models, defined in Sect. |2] and 
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Figure 1. Dark matter fraction /dm a.s a function of the shape parameter 
q for the SIS (lower symbols) and the Einasto (upper symbols) DM halo 
models, for two families with o"e8 = 150 and 300 km s~ for the spherical 
progenitors: the FO-built sub-families are in the top panel, while the EO- 
built ones are in bottom panel. The yellow, red, and black colours refer to the 
EO, E4, E7 model galaxies, respectively; symbols are filled for k = 0, and 
empty for k = \. The shape of the symbols (round or elliptical) indicates the 
FO or the EO view of a given model. See Sect. |3.3| and |3.4| for more details. 
The arrows indicate the trends of change of /dm for increasing flattening, 
and separate the symbols of the two DM profiles. 



summarized in Table |2] together with the main parameters charac- 
terizing the models. We next compare these temperatures with the 
observed X-ray properties of a sample of ETGs, extending to flat 
and rotating models the analysis carried out by Pll. In this Sec- 
tion, we take into account also the effect of a, that parametrizes the 
degree of thermalisation of the ordered motions. 

4.1 The effects of shape and stellar streaming motions on the 
model temperatures 

We explore here how T,, Ti^j and T^ depend on {q, k, a), i.e., galaxy 
flattening, rotational support, and degree of thermalisation of or- 
dered rotation. Three values of q = (1, 0.6, 0.3) are considered, that 
cover ETG morphologies from the spherical (EO) to the flattest ones 
(E7). The choice of the intermediate value q = 0.6 (corresponding 
to an E4) is motivated by the majority of ETGs having 0.55 < q < I 
(see Fig. |2j. Thus, the E7 models correspond to rare objects and 
represent quite an extreme behaviour, whereas the most common 
ETGs correspond to models with q between 1 and 0.6. We also 
consider the two extreme values of k: fully velocity dispersion sup- 
ported systems (k = 0), and isotropic rotators (k = 1); and three 
values of Q- = (0, 0.5, 1), in which respectively the pre-existing ISM 
has a null rotational velocity (all the stellar kinetic energy, includ- 
ing that of rotational motions, is thermalised, yth = l), or rotates 
with half the velocity of the stars (then 7,1, = 0.25), or has the same 
velocity as the stars (then no ordered stellar kinetic energy is ther- 
malised, 7th = 0, see Sect.|2ll. 
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Figure 2. Central stellar velocity dispersion as a function of the shape pa- 
rameter q for a sample of a 1400 nearby (z < 0.05) ETGs drawn from the 
SDSS DR4 (data taken form lNair & Abrahamll2010l) . 



Figure |3] shows T, for various descendants of spherical pro- 
genitors with cr^g = 150 and 300 Ism s"' (yellow circles), for two 
different DM haloes (SIS and Einasto in the left and right panels, 
respectively). As anticipated in Sect. 13.41 a major effect of flatten- 
ing is the decrease of cr^g of the descendants (with respect to the 
progenitor), that are then displaced on the left of their respective 
progenitor, in a way proportional to the flattening level |3 (arrows in 
Fig- El see also Figs.[3]and|4l(. 

The trend of T, with a pure change of shape in fully velocity 
dispersion supported models (k = 0), is due to the specific flatten- 
ing procedure (see Sect. |3.4t . In the FO-built sub-families (Fig. [3] 
top panels), flatter models are more concentrated than rounder ones, 
while in the EO-built sub-families (bottom panels), they are more 
extended and diluted. Thus pure flattening produces a different ef- 
fect on T,: in the FO-built cases T, increases (as T„; see below), 
whereas in the EO-built cases T, decreases. Therefore, we conclude 
that real flat galaxies can be either more or less bound than spherical 
galaxies of the same mass, depending on their mass concentration. 
Overall, however, the variation in T, for both sub-families is not 
large: the maximum variation, from the progenitor to the E7 model, 
is an increase of ~ 19 percent for the FO-built cases, and a decrease 
of ~ 18 percent for the EO-built ones. 

A larger effect on T, can instead be due to the presence of 
significant rotational support (empty symbols in Fig.|3}, if not ther- 
malised. In fact, when k = I, but a = (7th =1), the whole stellar 
kinetic energy, including the streaming one, is thermalised, and the 
T, values are coincident with those of the non-rotating case (full 
symbols), for the same galaxy shape. In the other cases of a i: 0, 
the rotational support always acts in the sense of reducing T„ and 
the flatter the shape, the larger can be the reduction. The strongest 
reduction of T, is obtained for an isotropic rotator (k = 1) E7 model, 
if the gas ejected from stars retains the same stellar streaming mo- 
tion (a = 1, 7th = 0): for the FO-built case, T, drops by ~ 50 
percent with respect to the EO model, and by ~ 60 percent with 
respect to the same E7 model with the ordered streaming motions 
fully thermalised (a = 0). For the EO-built case, T, drops by ~ 70 

^ In some works, instead of CTeg, the observations are used to measure the 
quantity Vntis = a/o"p + ^p. averaged within a central aperture (i.e., Re/8) 
by weighting with the surface brightness (see Appendix |bJ. For a chosen 
shape of the stellar distribution, Vrms = crp if i = 0, and whenever the 
galaxy is seen face-on. For any view, it can be shown that, for axisymmet- 
ric stellar distributions where the Satoh A:-decomposition is adopted, Vnns 
is independent of k; thus y,ms, c8 has a difi"erent behavior than (Teg , that is 
slightly lower iov k = 1 than for k = 0, for the edge-on view (Figs.[3]and|4). 
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Table 2. Summary of all parameters and temperatures. 



Symbol Meaning 



Intrinsic axial ratio of the stellar distribution: 0.3 ^ g ^ 1 

Satoh parameter, controls the amount of galaxy rotation: < /: < 1 

Degree of thermalisation of the ordered stellar motions v [Eq. (5)] 

Scahng factor between the ISM velocity and the stellar streaming motions (« = av), with < o- < 1 ; yth = ('!" ~ 1) 

Temperature equivalent of the thermalisation of the stellar motions and of the kinetic energy of SNe la events, for the unit mass of 

injected gas [Eq. (TJl 

Contribution to T^^j due to SNIa events [Eq. (6)]; it is regulated by a factor ;; < 1 (;; = 0.85 is generally adopted) 



1 
k 

7th 



TsN 
T, 
To- 



Contribution to T^^j due to stellar motions, defined by Eq. (2): T, = T^- + yth^'ro 

Contribution to T, due to stellar random motions, defined by Eq. ^3) 

Contribution to T, due to stellar ordered motions defined by Eq. (|4) 

Temperature equivalent of the change in gravitational energy of the injected gas, when flowing to the galactic centre [Eq. (8)] 

Temperature equivalent of the energy required to extract the unit mass of injected gas from the galaxy [Eq. l|9)]. If the injected gas 

rotates, then T^(a) = T^ - a^T^^t 

Mass averaged, subsonic escape temperature for the injected gas [Eq. UlU . It represents a fiducial upper limit to the observed 

temperature of outflows; T^"^ = 3/5 T^ 
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Figure 3. Distribution of T, for models derived from two spherical progenitors (yellow circles) of iTeg = 150 or 300 km s" , with SIS (left panel) or Einasto 
DM halo (right panel). Colours refer to the intrinsic flattening of the descendants: E4 (red) and E7 (black). At given intrinsic flattening (i.e., at fixed colour), 
the shape of the symbols indicates the l.o.s. inclination: EO (circle) or EO view (ellipse). Filled symbols refer to k = 0, empty symbols to k = 1. Three values 
of Q- are considered (a = 0,0.5, 1), as indicated for the FO-built sub-families of the iTe8=300 km s" progenitor When k = 1, the increase of a decreases T, at 
any a-^g, from the value coincident with the non-rotating case (k = 0) when a = 0, down to values that are lower for larger flattenings, when a = 0.5 and 1. In 
each sub-family, models with same intiinsic shape have the same T,, independent of the EO or EO view, for fixed (k, a). 



percent with respect to the EO model, and by ~ 60 percent with 
respect to the same E7 model with a = 0. These percentages are 
obviously extreme values; the T, reduction is lower for milder flat- 
tenings, and for k and a values smaller than 1. For a fixed galaxy 
shape, all possible (k, a) combinations fill a sort of triangular area 
on the (cTeg, T,) plane, identifiable by linking the symbols of a given 
q (colour). Clearly, the rounder the galaxies, the weaker the effect 
of k and, consequently, of a variations. All the above effects are in- 
dependent of the galaxy luminosity (mass), and both the a^g = 300 
and 150 km s"' families show the same (rescaled) behaviour in the 
(cre8,r.) plane. 

The trends described above are independent of the specific 
DM halo profile: models with an Einasto DM halo (right panel) 



show the same pattern as those with a SIS DM halo (left panel), 
just with a different normalisation due to the larger total DM con- 
tent (see Table [TJ. Similar results hold also for the Hemquist and 
NEW DM haloes. 

In Eig.|4]we plot T^ for the same families in Eig.|3] Similarly 
to what happens for T,, for fully velocity dispersion supported mod- 
els, T„ gets larger with flattening for the EO-built sub-family (~ 21 
percent), while it decreases for the more diluted EO-built models 
(~ 12 percent). Stellar streaming, when a > 0, acts in the sense 
of making the gas less bound, due to the centrifugal support of the 
injected gas, and then T^ decreases with increasing a, at any fixed 
flat shape. This effect is maximum when a = I and the gas rotates 
as the stars. However, the decrease in T^ due to galaxy rotation is 
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Figure 4. T^ for the same models in Fig. |3](SIS and Einasto DM halo in the left and right panels, respectively); the same notation appUes. The shaded area 
shows the range spanned by T[„j, considering the variations in T, in Fig.|3] for ;; = 0.85. 



lower than obtained for T,: for both sub-families, T^ drops at most 
by ~ 13 percent (for the E7 models), between the two extreme 
cases of Q' = and a = I. This produces that, in the FO-built case, 
T„ keeps always larger than for the progenitor when q decreases, 
even for k = a = 1, while T, of rotating galaxies could become 
significantly lower than for the progenitor. Note that there are two 
compensating effects from stellar streaming when a t 0: the stel- 
lar heating is lower than for k = 0, but the gas is also less bound. 
Also these results are independent of the DM halo profile, as can be 
judged from the right panel of Fig.|4] that refers to the same models 
of the right panel of Fig.|3] 

Since the flatter is the galaxy, the more it can be rotationally 
supported (and the more is rotationally supported, the larger is the 
effect of a corotating ISM), the effect of rotation is dependent on 
the degree of flattening, and thus it may prove difficult to disentan- 
gle observationally the two distinct effects due to shape and kine- 
matics. On the theory side, we recall that a simplifying assumption 
made here is that u = av, and yth = (« - 1)^, while in reality the 
kinematical difference between stars and pre-existing gas may be 
more complex, as the extent of thermalisation; only numerical sim- 
ulations will be able to establish what are t he net effects on th e gas 
evolution of stellar streaming motions (see lNegri et al. 11201 3 al ibi). 

Figure |4] finally shows the well known fact that the con- 
tribution from SNIa's dominates the gas injection energy, since 
Tinj >> T,. This contribution (i.e., Tsn) is independent of galaxy 
mass, which results into lower-mass models having T^ far lower 
than Tinj, and T^ reaching Tinj for 200 km s"' < cr^g < 250 km s"', 
depending on the DM profile. Galaxies with (T^g < 200 km s"' con- 
sequently are more prone to an outflow, and then to have a low hot 
gas content, as already suggested in the past by numerical s imula- 
tions and by observations (CDPR, Sarazin e t al. 2001 ; Davi d et alj 
bOOeLIPellegrini et al.ll2007l;lTrinchieri et al.l boOS). These findings 
are based on the assumption of a high thermalisation efficiency for 
SNIa's (rj = 0.85), and the quoted cr^^ c ritical values become lower 
for lower rj values (as indicated by, e.g.. lThomton et al.lll998l) . that 



decrease Ti„j. Variations in the dark matter may alter the T„ values, 
but small changes in T^ are found here for differences in the DM 
profile, and possible variations in the total DM amount cannot be 
very large, given the constraints from dynamical modelings within 
R^, and from cosmological simulations (taken into account here. 
Sect. 13.3b . Indeed, the results for the Hemquist families are essen- 
tially identical to what we have shown for the Einasto halo, whereas 
for the NFW profile the trends are the same, but all the temperatures 
are shifted to higher values, due to their larger amounts of DM. 



Finally we comment on the preliminary investigation about 
the role of flattening and rotation on the global energetics of 
the ISM in CP96. They built fully analytical axisymmetric two- 
component galaxy models, where the stellar and dark mass distri- 
butions were described by the Miyamoto-Nagai potential-density 
pair. They varied the shape of both the stellar component and the 
dark matter halo from flat to spherical, and the amount of azimuthal 
ordered motions through the Satoh fe-decomposition. Their conclu- 
sion was that, for quite round systems, flattening can have a sub- 
stantial effect in reducing the binding energy of the hot gas, con- 
trary to galaxy rotation that seemed to have a negligible role. The 
opposite was suggested for very flat systems. These results were 
confi rmed by 2D hydrodynamical simulations dD'Ercole & Ciottil 
1 19981) . We stress here that the models in CP96, while capturing 
the main effects of flattening and rotation on the global energet- 
ics of the hot haloes of ETGs, were not tailored to reproduce 
in detail the observed properties of real ETGs. Our current find- 
ings, based on more realistic galaxy models, reveal a more compli- 
cated situation, where the flattening importance on the ISM status 
is mediated by the amount of rotation and its specific theiTnalisa- 
tion history. Remarkably enough, however, when flattening a two- 
component Miyamoto-Nagai model following the same procedure 
here adopted (constant M,, R^ and Mh), its T^ remains close to that 
of its spherical progenitor (i.e., the model moves almost parallel to 
the solid lines in fig. 2 in CP96). 
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Table 3. Observed Properties of the ETG Sample with X-ray properties for the hot gas from Chandra observations. 



Name 


d 


log(LK) 


k^Tx 


Lx 


•'max 


O-eS 


Vmax/o"c8 


Ref 


Type 


1 




(Mpc) 


(iKo) 


(keV) 


(lO'tOergs-l) 


(kms"') 


(kms"') 






RC3 


2MASS 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



NGC 720 


27.6 


11.31 


0.54 


5.06 


100 


241 


0.41 


Binney et al. 1990 


E5 


0.55 


NGC 821 


24.1 


10.94 


0.15 


2.13x10^3 


120 


200 


0.60 


Coccato et al. 2009 


E6? 


0.62 


NGC1023 


11.4 


10.95 


0.32 


6.25x10^2 


250 


204 


1.23 


Noordemieer et al. 2008 


SBO 


0.38 


NGC1052 


19.4 


10.93 


0.34 


4.37x10-1 


120 


215 


0.56 


Milone et al. 2008 


E4 


0.70 


NGC1316 


21.4 


11.76 


0.60 


5.35 


150 


230 


0.65 


Bedregal et al. 2006 


SABO 


0.72 


NGC1427 


23.5 


10.82 


0.38 


5.94x10^2 


45 


171 


0.26 


D'Onofrio et al. 1995 


cD 


- 


NGC 1549 


19.6 


11.20 


0.35 


3.08x10"' 


40 


210 


0.19 


Longo et al. 1994 


EO-1 


0.90 


NGC2434 


21.5 


10.84 


0.52 


7.56x10-' 


20 


205 


0.10 


Carollo et al. 1994 


EO-1 


0.98 


NGC2768 


22.3 


11.23 


0.34 


1.26 


195 


205 


0.95 


Proctor et al. 2009 


E6 


0.46 


NGC3115 


9.6 


10.94 


0.44 


2.51x10-2 


260 


239 


1.09 


Fisher 1997 


SO 


0.39 


NGC3377 


11.2 


10.45 


0.22 


1.17x10-- 


97 


144 


0.67 


Simien et al. 2002 


E5-6 


0.58 


NGC3379 


10.5 


10.87 


0.25 


4.69x10-2 


60 


216 


0.28 


Weijmans et al. 2009 


El 


0.85 


NGC3384 


11.5 


10.75 


0.25 


3.50x10-2 


150 


161 


0.93 


Fisher 1997 


SBO 


0.51 


NGC3585 


20.0 


11.25 


0.36 


1.47x10-' 


200 


198 


1.01 


Fisher 1997 


E6 


0.63 


NGC3923 


22.9 


11.45 


0.45 


4.41 


31 


250 


0.12 


Norris et al. 2008 


E4-5 


0.64 


NGC4125 


23.8 


11.35 


0.41 


3.18 


150 


227 


0.66 


Pu et al. 2010 


E6 pec 


0.63 


NGC4261 


31.6 


11.43 


0.66 


7.02 


50 


300 


0.17 


Bender et al. 1994 


E2-3 


0.86 


NGC4278 


16.0 


10.87 


0.32 


2.63x10-' 


60 


252 


0.24 


Bender etal. 1994 


El-2 


0.93 


NGC4365 


20.4 


11.30 


0.44 


5.12x10-' 


80 


245 


0.33 


Surma etal. 1995 


E3 


0.74 


NGC4374 


18.3 


11.37 


0.63 


5.95 


60 


292 


0.21 


Coccato et al. 2009 


El 


0.92 


NGC4382 


18.4 


11.41 


0.40 


1.19 


70 


187 


0.37 


Fisher 1997 


SAO 


0.67 


NGC4472 


16.2 


11.60 


0.80 


18.9 


73 


294 


0.25 


Fisher etal. 1995 


E2 


0.81 


NGC4473 


15.7 


10.86 


0.35 


1.85x10-' 


70 


192 


0.36 


Emsellem et al. 2004 


E5 


0.54 


NGC4526 


16.9 


11.20 


0.33 


3.28x10-' 


246 


232 


1.06 


Pellegrini et al. 1997 


SABO 


0.43 


NGC4552 


15.3 


11.01 


0.52 


2.31 


17 


268 


0.06 


Krajnovic et al. 2008 


EO-1 


0.94 


NGC4621 


18.2 


11.16 


0.27 


6.08x10-' 


140 


225 


0.62 


Bender et al. 1994 


E5 


0.65 


NGC4649 


16.8 


11.49 


0.77 


11.7 


120 


315 


0.38 


Pinkney et al. 2003 


E2 


0.81 


NGC4697 


11.7 


10.92 


0.33 


1.91x10-' 


115 


174 


0.66 


De Lorenzi et al. 2008 


E6 


0.63 


NGC5866 


15.3 


10.95 


0.35 


2.42x10-' 


210 


159 


1.32 


Neistein et al. 1999 


SAO 


0.42 



Notes: (1) Galaxy name. (2) Distances from BKF ITonrv et al.l200lh . (3) Logarithm of the K band luminosity, assuming Kq = 3.33 mag, taken from BKF. 
(4) - (5) Hot gas temperature and the 0.3 - 8 keV gas luminosity, from BKF. (6) - (8): maximum velocity of rotation, stellar velocity dispersion, as the 
luminosity-weighted average within an aperture of radius Re/8 (from PI 1), and their ratio. (9) References for Vmax in column (6). (10) Morphological type 
from RC3. (11) Axial ratio in the K, band, from 2MASS. 



4.2 Comparison with observed ETGs properties in the 
X-rays 

Now we compare our estimates for 7", and T^ with observed tem- 
peratures Tx and gas content, respectively, for the BKF sample. Ob- 
served X-ray and A'-band luminosities, central stellar velocity dis- 
persions, rotation velocities, and galaxy shape are given in Tab. [3] 
Figure |5] shows observed (points) and model temperatures (lines) 
versus Ly^, with Tx plotted with a different colour reflecting the 
ETG shape and rotational support. In order to make the compari- 
son between models and observations more consistent, the Lk of 
the models have been calculated using the mean Lk/Lv = 3.4 of 
the ETGs in this sample. When comparing observed and model 
temperature values at fixed Lk, the mass is the same for all the 
models, and it should be roughly so also for the observed ETGs. At 
fixed Lk, then, the r» variation due to shape and stellar streaming is 
obtained from all the T, values of the descendants of a progenitor 
(we consider the Einasto models of Fig. |3}, regardless of the FO 
or EO view. When comparing observed and model temperatures as 
a function of a-^^ instead, at any CTeg the mass could be different. 
Figure |5] shows again how the effect of rotation can be potentially 
stronger than that of shape (as already indicated by Fig.|3): moder- 
ately larger or smaller T, can be produced by flattening, depending 
on the way it is realized, while T, can be much lowered by galactic 
rotation. Thus, at fixed Lk, the largest variation in T, with respect 



to the spherical case (yellow line) does not come from a variation 
of shape, but is a decrease of T, due to rotation with a = \ (no 
thermalisation). 

In Fig. |5] all Tx are larger than T,, yet much closer to T, than 
to Linj ~ 1.5-2 keV (see Fig.O. This may be evidence of two 
facts: either the SNIa's thermalisation rj is low, or the gas flows es- 
tablish themselves at a temperature close to the virial temperaturqj 
of the galaxy, and most of the SNIa's input is spent in cooling (in 
gas-rich ETGs), or in lifting the gas from the potential well, and 
imparting bulk velocity to the outflowing gas (in gas-poorer E TGs; 
see Pll for a q uantification of these effects, and Tang et al. 20091 
and lLi et al.ll2011al for addressing also other solutions to this prob- 
lem). A combination of the two explanations may also be at place, 
of course. In any case, the proximity of the observed Tx values to 
r, provides an empirical evidence of the importance of the study 
of T, and its variations; it would not have been so, if we had found 
Tx to be closer to ri„j. 

Going now into the question of whether possible effects from 
shape and stellar streaming are apparent on Tx, Fig. [5] shows a 



^ The integrals in the definitions of To- and T^at (Eqs. 3 and 4) are also 
used to compute the total kinetic energy of the stellar motions that enters 
the virial theorem for the stellar component; thus the mass weighted tem- 
perature r», with 7ii, = 1, is often referred to as the gas virial temperature. 
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Figure 5. T, (lines) for all the progenitors in Tab. [T]with the Einasto halo, and their descendants (with a = 0ov a = 1), and the observed Tx (circles), for the 
ETGs in the BKF sample, as a function of Lk. The yellow line refers to the progenitor ETGs, the red lines to the E4 shape, the black ones to the E7 shape; 
lines are soUd for k = 0, and dashed for k = \ and a = I. For each colour, lines representing the case of k = \ and ff = are coincident with the solid lines. 
Thick and thin lines refer to the FO and EO-built models, respectively. Left panel: the colour-coding for the observed ETGs indicates their ellipticity e = \- q, 
as measured in the Ks-band, from 2MASS (see Tab. [5), and is cahbrated as for the previous figures [i.e., to be yellow for the EO (c/ = 1, e = 0), red for the 
E4 (q = 0.6, e = 0.4), and black for the E7 {q = 0.3, e = 0.7); see the colour bar on the right for the e of the other colours]. Right panel: the colour-coding 
indicates the rotational support Vmax/o"e8. Note that, for an EO view, the models of the dashed lines would have Vmax/o"c8 = 0.90 - 0.95 for the E4 case, and 
^max/o"e8 = 1.5 - 1.6 for the E7 case. 



mild indication that flatter shapes and more rotationally supported 
ETGs tend to show a lower Tx, with respect to rounder, less rotating 
ETGs. This is similar to the recent result by S13, that fast rotators 
seem to be confined to lower temperatures than slow rotators. S13 
suggested that in ETGs with a larger degree of rotational support, 
the kinetic energy associated with the stellar ordered motions may 
be thermalised less efficiently. Preliminary results of h ydrodynam- 
ical si mulations seem to indicate that this is the case JNegri et alj 
l2013al lbh. Based on the analysis of Sect. 14.11 we can propose an 
explanation for the trends in Fig.|5]that seems consistent with the 
present data, and that represents a prediction for when the sample 
in Fig.[5] still quite small, will be hopefully enlarged: stellar stream- 
ing always reduces r,, in a way proportional to k and a, up to an 
amount that can be as large as 60 percent, while flattening itself is 
less important (though necessary for rotation to be important). 

We now move to consider possible effects from shape and rota- 
tion on the gas content, starting from the hypothesis that it is linked 
to the relative size of T^ and ri„j. Figure |6] shows the observed Lx 
values, together with the energies (referring to the mass of gas in- 
jected in the unit time) describing the stellar heating L,, the SNIa's 
heating Lsn> and the requirement for escape L^. All quantities are 
normalised to Lk, and are computed using the expression for the 
stellar mass loss rate M, and the temperatures defined in Sect. |2] 
such that L- = M.R- = 3kj^M,T-/(2i^mp), L, = 2,k^M,T.I(2iimp), 
and LsN = J/L^sn^sn- Figure |6] clearly shows that flatter and more 
rotationally supported ETGs tend to have a lower Lx/Lk, as already 
known (see Sect. [T). As shown in Sect. 14. 1 1 however, the effect of 
shape or rotation on L^ is sinall, so we cannot claim an important 
direct role for these two major galactic properties on determining 
a lower gas content and then Lx. It has been suggested that a pos- 
sible indirect effect could come from galactic rotation if it is effec- 
tive in creating a gas disc, where gas cooling is triggered, the tem- 



perature is lowered, and then Lx is reduced (IBrighenti & Mathewsl 
119971) . Other possibilities may be related to global instabilities of 
rotating flows, perh aps associated with inefficient thermalisation 
JNegri et al.ll2013bl) . Note that the X-ray emissivity is also depen- 
dent on the gas temperature, and one could think that the lower Lx 
of flat/rotating ETGs could be due to having these preferentially a 
lower Lx (Fig.[5]l. This cannot be the explanation, though, because 
the emissivity in the 0.3-8 keV band decreases very mildly with de- 
creasing temperature, for temperatures below a value of ~ 1 keV. 
Another explanation for a lower Lx/Lk could be a lower stellar age 
in fast rotators: indeed, S13 found that molecular gas and young 
stellar populations are detected only in fast rotators across the en- 
tire ATLAS^'^ sa mple, and a younger ag e is known to be linked to 
a lower Lx (BKE lO'Sullivan et allboOli) . 

Also, note in Fig. [6] how there are many ETGs with Lx lower 
than L,: they do not even radiate L,. For them, the outflow must be 
very important, and must have employed almost all of Lsn. Numer- 
ical simulations (CDPR) have already shown that Lx can be even 
lower than L, during winds/outflows. As L^ increases and becomes 
closer to Lsn, the ETGs below L, disappear. 

Finally, BKF found a positive correlation between Lx and Lx, 
valid down to the gas-poor galaxies, that on average have the shal- 
lowest potentials (Lk ^ 2x 10"Lko and o"e8 5 200 km s"') and are 
expected to host outflows; also in the Lx - Lk plot, in the large 
variation in Lx at any Lk <; 7 x 10"'Lko5 ETGs with the low- 
est Lx tend to be the coldest ones. This seemed a puzzle, if out- 
flows (lower Lx) are expected to be hotter than inflows (large Lx). 
Pll re-examined the Lx behavior, considering Lx values rescaled 
by r,, where for the latter the case of a spherical galaxy with an 
isotropic velocity dispersion tensor was taken. The gas of ETGs 
with cr^f, < 200 km s"' turned out to be colder in an absolute sense, 
but to have the largest values for the rescaled Lx; the latter also 
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Figure 6. Comparison between the ran of various luminosities vs. Lk, for all the progenitors in Tab.[T](Einasto halo case) and a selection of their descendants 
(only cases of ff = 0, 1), and for the observed Lx for the ETGs in the BKF sample; all quantities are normalised to Lk. The lower group of lines gives L,/Lk, 
the upper one gives Lg /Lk, the horizontal line is Lsn/^-k for r\ = 0.85. The yellow line refers to the progenitor ETGs, the red lines to the E4 shape, the black 
ones to the E7 shape; lines are solid for k = 0, and dashed for k = I and a = I, respectively. Thick and thin lines refer to the FO and EO-built models, 
respectively. The colour-coding for the observed ETGs indicates their ellipticity e = I- q (left panel) and the rotational support V,„ax/o"e8 (right panel), and is 
calibrated as for the previous Fig.|5] 
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Figure 7. Relative temperatures versus iXeg for obsei'ved ETGs and models (only EO-built models observed EO; the plot is very similar' for the FO-built 
models). The colours and linetype for the models are the same as in Fig.|5] T, values (lower bundle of lines) ai'e normalised to r» corresponding to a spherical 
galaxy with an isotropic velocity tensor (given by the yellow line in Fig.|5); 7"™^!" (upper bundle of lines) and Tx's are normalised to the T, of the corresponding 
shape (neglecting rotation). In particular, observed ETGs with q in the range [1,0.75] are normalised to T, of the EO model, those with q in the range [0.75,0.45] 
are normalised to T, of the E4 model, whereas those with q < 0.45 are normalised to T, of the E7 model. The colour-coding for the obsei'ved ETGs indicates 
their logLx(erg s" ), in the left panel, and the rotational support Vmax/o"c8, in the right panel, as indicated by the colour bar on the right. 



tend to the (rescaled) T™^!", a fiducial average temperature of gas 
in outflow (Sect. I2.2I >. consistent with the expectation for the gas 
flow status of these ETGs. At intermediate cr^g values, though (200 
km s" < iTe8 < 250 km s" ), still the rescaled Tx seemed to be 
lower for the gas poorest ETGs, which remained unexplained. Can 
we suggest that at these intermediate cr^s the effect of shape or ro- 



tation is responsible for a decrease of both Tx and Lx? And then, 
when normalising Tx by a T, appropriate for the shape and rota- 
tional level of the host galaxy, the result is larger than when nor- 
malising by the T, of a spherical ETG of same cr^gl Figure|7]shows 
Tx values normalised by a T, appropriate for the apparent shape of 
the galaxy, neglecting rotation; it confirms that at low cr^g the gas 
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is relatively hotter, and reaches close to T^'^jT, (Pll). In the right 
panel, the colour indicates the level of rotational support Vmax/Ces; 
if rotation is not neglected, in this plot the properly normalised Tx 
of rotating ETGs would become higher (see, e.g., Fig.|5}- Indeed, 
in the intermediate cr^% region, some of the gas-poor ETGs with the 
lowest Tx/T,, are also highly rotating, and then their position in 
the plot would be higher, thus at least partially accounting for the 
segregation. 



5 SUMMARY AND CONCLUSIONS 

In this paper we investigated the relationship between the tempera- 
ture and luminosity of the hot X-ray emitting haloes of ETGs and 
the galactic shape and rotational s upport. This work is an extension 
of previous similar studies (CP96. |Pellegrini et al.lll997l Pll). 

By solving the Jeans equations, we built a large set of axisym- 
metric three-component (stars, dark matter halo, SMBH) galaxy 
models, representative of observed ETGs. We varied the degree of 
flattening and rotational support of the stellar component, in order 
to establish what is the dependence of the temperature and binding 
energy of the injected gas on the observed galaxy shape and internal 
kinematics. For the injected gas, we defined the equivalent tempera- 
ture of stellar motions as T, = T„- + ythTiot! where the parameter yth 
takes into account how much of the ordered rotation of the galaxy 
is eventually thermalised by the stellar mass losses. We considered 
the simplified case in which the pre-existing gas velocity is propor- 
tional to the stellar streaming velocity. When pre-existing gas and 
stars rotate with the same velocity, no ordered stellar kinetic en- 
ergy is thermalised (7,1, = 0); when the gas is at rest, all the kinetic 
energy due to stellar streaming is thermalised (y,h = 1). We also 
defined the temperature equivalent (T„) to the binding energy for 
the injected gas. Our main results are as follows: 

• The major effect of flattening the stellar component of a spher- 
ical ETG is a decrease of the observed a-^s value of its flatter coun- 
terparts of same mass and same circularized R^. This decrease is 
proportional to the flattening level and depends also on the view- 
ing angle. For each shape, in velocity dispersion supported models, 
the decrease is larger for the FO-view and smaller for the EO-view, 
and it reaches ~ 35 percent for the E7 shape seen FO. In isotropic 
rotators, the decrease is instead larger for the EO view. Thus, ETGs 
with the same roundish appearance and the same observed a^s may 
have a significantly difl'erent mass. This finding raises the issue of 
the reliability of the use of a-^s as a proxy for the dynamical mass 
of an ETG, a point particularly relevant for studies of the hot haloes 
properties, that mainly depend on the galaxy mass. 

• Flatter models can be either more or less concentrated than 
rounder ones of the same mass and same circularized R^, depend- 
ing on how they are built: if R^ is kept constant for a FO view, flatter 
models are more concentrated and bounded than the round counter- 
part; the opposite is true if R^ is kept constant for an EO view. As a 
consequence, the effect of a pure change of shape is an increase of 
r, and r^ in the first way of flattening, and a decrease in the sec- 
ond one. Overall, however, the variation in T, for both cases is mild, 
within ~ 20 percent even for the maximum degree of flattening (the 
E7 model). Similarly, T^ gets larger by at most ~ 20 percent, and 
decreases by at most ~ 12 percent. 

• A more significant effect on T, can be due to the amount of 
rotational support. The isotropic rotator case is investigated here 
(k = I). If 7(h = 1, the whole stellar kinetic energy, including the 
streaming one, is thermalised, and T, coincides with that of the 



non-rotating case, for the same galaxy shape. If yn, < 1, T, is in- 
stead always reduced, and the larger so the flatter is the shape. The 
strongest reduction is obtained for the E7 models when 7,11 = 0, and 
then T, can drop by 50 - 70 percent with respect to the EO models. 
Thus, the presence of stellar streaming, when not thermalised, acts 
always in the sense of decreasing T„ and the size of this decrease 
depends on the relative motion between pre-existing gas and stars. 
Clearly, the flatter the galaxies, the stronger can be the rotational 
support, and, consequently its potential effect on T,. Thus the ef- 
fect of rotation is dependent on the degree of flattening, and, as a 
minimum, it requires a flat shape as a premise. 

• Since stellar streaming acts in the sense of making the gas less 
bound due to the centrifugal support, at any fixed galaxy shape and 
rotational support T„ decreases in proportion to how the velocity 
field of the ISM is close to that of the stars. However, this decrease 
is lower than that produced on T,: T„ drops at most by ~ 13 percent 
(for the E7 models), between the two extreme cases of ISM at rest 
and gas rotating as the stars. Note then two compensating effects 
from stellar streaming when yth < 1: T. is lower than for a non- 
rotating galaxy, but the gas is also less bound. 

• All the above trends and effects are independent of the galaxy 
luminosity (mass), and dark halo shape. Only the normalization of 
T, and T^ changes, if the dark halo mass changes (e.g., it increases 
from the SIS to the Einasto to the Hemquist to the NEW halo mod- 
els). 

The comparison of the above results with observed Tx and Lx 
for the ETGs in the Chandra sample of BKF shows that: 

• All observed Tx are larger than T„ but much closer to T, than 
to TsN- T, ranges between 0.1 and 0.4 keV, for 150 km s"' < iXeg < 
300 km s"' (the lower end being possibly even lower, depending 
on the effects of rotation), while Tsn ~ 1.5 keV (for a thermali- 
sation parameter r) = 0.85), so that the contribution from SNIa's 
dominates the gas injection energy (Tjnj ~ 1.5-2 keV), that is then 
practically insensitive to changes in T, due to the galaxy shape or 
kinematics. The proximity of the observed Tx values to T, indicates 
that rj may be lower, and/or that the gas tends to establish itself at 
a temperature close to the virial temperature, in all flow phases. In 
any case, this proximity provides an empirical confirmation of the 
relevance of a study of r» and its variations. 

• For 200 km s"' < cr^g < 250 km s"', T„ becomes larger than 
Tjnj, and inflows in these galaxies can become important. These iXeg 
values could be lower if 77 < 0.85, and if the dark matter amount is 
larger than assumed here. Galaxies with a^^ < 200 km s"' (Lk ^ 
2 X 10"Lko), instead have a high probability of hosting an outflow, 
and then a low hot gas content, as confirmed by observations of 
LxZ-i-K vs. Lk. The ratio Lx/Lk seems to be lower also for a flatter 
shape and larger rotation. However, the effect of shape or rotation 
on T„ is small, thus these two major galactic properties are not 
expected to play an important direct role in determining the gas 
content. Indirect effects may be more likely (as galactic rotation 
triggering large-scale instabilities in the gas, or the lower age of 
fast rotators). For many ETGs, Lx/^-k is much lower than L,/Lk, 
indicating clearly the presence of an outflow. 

• We find a mild indication that, at fixed Lk, flatter shapes and 
more rotationally supported ETGs show a lower Tx, with respect 
to rounder, less rotating ETGs (similarly to what recently found by 
SI3). This tendency can be easily explained by the effects predicted 
here on T, due to flattening and rotation. Since, for a fixed galaxy 
mass, a decrease of T, due to rotation is predicted to be potentially 
stronger than produced by shape without rotation, we propose that 
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not thermalised stellar streaming is a more efficient cause of the 
(possibly) lower Tx- 

• Extending a Pll result, we find that, when rescaled by a 
T, value representative of a fully velocity dispersion supported 
ETG of the corresponding shape, the TxIT, values of outflows (at 
CTeg < 200 km s"') are larger than those of inflows (at iXeg <; 250 
km s"'). At 200 km s"' < cr^^ < 250 km s"', where the T^IT, val- 
ues seem to be lower for the gas poorest ETGs, part of the trend 
could be explained by a few of these ETGs being highly rotating. 
Note that the observed galactic rotation depends on the viewing an- 
gle at which these low Tx and low Lx ETGs are seen, thus it is 
difficult to find a clear indication of a systematic trend in the data. 

As a necessary parallel investigation, we are going to perform 
2D hydrodynamical simulations to exp lore the actual be haviour of 
the ISM for the models here studied (iNegri et al.ll2013bh . The goal 
will be to establish, as a function of galactic structure and kine- 
matics, the presence of major hydrodynamical instabilities, the ef- 
ficiency of thermalisation of stellar streaming motions, and their 
impact on Tx and Lx. 
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APPENDIX A: THE FLUID EQUATIONS IN THE 
PRESENCE OF SOURCE TERMS 

The equations of fluid dynamics in the presence of different sources 
of mass, momentum and energy, under the si mplifying assumption 
of isotropy of the mass losses, can be written ( ID'Ercole et alJlOOOh 

as 
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Here ./^(x, t) = S^ is the total mass return per unit time and vol- 
ume, due to different sources (e.g., stellar winds and SNla events) 
associated with a stellar population with streaming velocity v(x, t) 
and velocity dispersion cr-(x, t). e,(x, t) is the internal energy return 
per unit mass and time of the i-th source field, and ihj(x, t) is the 
modulus of the relative velocity of the material injected by the i- 
th source field with respect to the source star. Finally, J^ are the 
bolometric radiative losses per unit time and volume. 

For example, in applications as the one in this paper, ^ is 
represented by the sum of stellar winds and SNla explosions ejecta. 
Often, s„i„d and Ms.wind are neglected, being significantly smaller 
than the contribution of Tr(a"^), while the opposite holds for SNla 
events, where ./i^sNTr(a"-) is negligible with respect to the energy 
injection due to the SNla explosions. In our case, yth in Eq. ^ 
derives from the full mass injection term in Eq. ( |A3) . where both 
contributions are taken into account. 



APPENDIX B: THE CODE 

All the relevant dynamical properties of th e models are com- 
puted using a numerical code built on purpose jPosacki et al.l2013l : 
|Posackj||2011l) . Starting from an axisymmetric density distribution 
the code computes in terms of elliptic integrals the associated grav- 
itational potential and vertical and radial forces. The code then 
solves the Jeans equations in cylindrical coordinates (R,z,ip)- For 
an axisymmetric density distribution p,(R,z) supported by a two- 
integral phase-space distribution function (DF), the Jeans equations 
write 
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dz 
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(e.g. lSinnev & Tremaindl 1987b . where 0,oi is the sum of the grav- 
itational potentials of all the components (e.g. stars, dark halo, 
black hole). As well known, for a two-integral DF (1) the veloc- 
ity dispersion tensor is diagonal and aligned with the coordinate 
system; (2) the radial and vertical velocity dispersions are equal, 
i.e. (Tji = a- = cr; (3) the only non-zero streaming motion is in the 
azimuthal direction. 



In order to control the amount of ordered azimuthal velocity 
v^, we adopted the ^-decomposition introduced bv lSatohl ( 119800 



k^iv 
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and then it follows 



0-2 + (1 - l^)(vl - a\ 
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where < A: < 1 . The case k = \ corresponds to the isotropic rota- 
tor, while for /: = no net rotation is present and all the flattening 
is due to the azimuthal velocity dispersion cr^. In principle, k can 
be a function of {R, z), and so more c omplicated (realisti c) velocity 
fields can be realized (CP96; see also lNegri et al.l2013ah . The code 
then projects all the relevant kinematical fields, together with the 
stellar density. 

The projections along a general l.o.s. of the stellar density p., 
streaming velocity v and velocity dispersion tensor cr^ are 
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respectively, where (,> is the scalar product, n is the l.o.s. direction 
and / is the integration path along n. Note that if a rotational sup- 
port is present, then cTp is not the l.o.s. (i.e. the observed) velocity 
dispersion a^^^, given by 



o-f + Vv - vL' 
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where Vp is the projection of {n, v) (CP96). In particular, the face- 
on projections are 
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where all the integrations are performed at fixed z- 

After the projection, the code calculates the corresponding cir- 
cularized effective radius R^. In practice, the stellar projected den- 
sity is integrated on the isodensity curves, and Rl = q\o!,a\, where 
fle and ^los are the semi-major axis and l.o.s. axial ratio of the el- 
lipse of half projected luminosity, respectively. In general , ^i^s is 
a function of the l.o.s. inclination ( ILanzoni & Ciottil 120031) : for a 
face-on projection ^los = 1, while for the edge-on case q\o^ - q (see 
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Sect. 13. lb . Then, the corresponding luminosity averaged aperture 
velocity dispersion 

is computed within a circular aperture of R^l%. 

For comparison with other works, we followed also the ap- 
proach of the ATLAS ^'^ project (Cappellari et al. 20 11 , .20 IZ) . cal- 
culating the quantity 

VL = ^? + Vp2 = ^L + vL (B16) 

and its corresponding luminosity averaged mean within i?e/8, ac- 
cording to a definition analogous to Eq. ( IB15t . Finally, for a given 
model, the code evaluates all the volume integrals in Sect. |2l by 
using a standard finite-difference scheme. 
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